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Abstract. In this paper, we obtain the analytical solutions of two kinds of 
transcendental equations with numerous applications in college physics by means 
of Lagrange inversion theorem, and rewrite them in the form of ratio of rational 
polynomials by second order Pade approximant afterwards from a practical and 
instructional perspective. Our method is illustrated in a pedagogical manner for the 
purpose that students at the undergraduate level will be beneficial. The approximate 
formulas introduced in the paper can be applied to abundant examples in physics 
textbooks, such as Fraunhofer single slit diffraction, Wien’s displacement law and 
Schrodinger equation with single or double 6 potential. These formulas, consequently, 
can reach considerable accuracies according to the numerical results, therefore they 
promise to act as valuable ingredients in the standard teaching curriculum. 
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1. Introduction 

Transcendental equations of certain types frequently emerge in seemingly unrelated 
branches of physics, and the roots of these equations appear in a great deal of 
applications. Furthermore, the ascertainment of zeros of equations is a problem 
commonly encountered in a broad spectrum of scientihc applications. A wide variety 
of root hnding algorithms [1] are available to approximate the solutions to any desired 
degree of accuracy though, exact analytical solutions to physical problems, which provide 
better insights into the physical signihcances of associated parameters than purely 
numerical solutions, are always desirable and preferable. 

The approximation formula, however, plays an unique role in teaching physics since 
it hlls the gap between analytical approach and numerical solution. Frustratedly, many 
useful approximation techniques, aiming to solve special physical problems, are not 
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really familiar to college students. Quite recently, Kevin Rapedius presented complex 
resonance states (or Siegert states) that describe the tunnelling decay of a trapped 
quantum particle by means of Siegert approximation method [2]. Augusto Belendez et al 
obtained a simple but highly accurate approximate expression for the period of a simple 
pendulum at the aid of Kidd-Fogg approximate formula [3]. In this paper, we intend 
to draw attention to the Lagrange inversion theorem [4,5] and Fade approximant [6] by 
two kinds of transcendental equations that have made their appearances in a variety of 
applications. The hrst kind of equation is 

tanx = Kx, (la) 

cot X = KX, {lb) 

which, splendidly, can be used to hnd the correlation to the spring oscillations due to the 
non-ignorable mass of spring [7], to determine the positions of maxima of Fraunhofer 
single slit diffraction [8] and the eigenvalues of inhnite square potential well with a 
residual 5-function interaction [9,10]. The second one is 

W{x)e^^^'^ = X (2) 

where W{x) is the celebrated Lambert W function [11], which has undergone a 
renewed interest during the past three decades and owns a diversity of multidisciplinary 
applications. Numerous quantities, such as Wien’s displacement constant [12], can 
be expressed in closed form in terms of the Lambert W function. Here we refer the 
motivated readers to [13,14] and references therein for more information. 

The main purpose of this paper is to provide a pedagogical derivation of the 
analytical solutions of transcendental equations (1) and (2) with the help of Lagrange 
inversion theorem, and bother Fade approximant subsequently to express them in 
element forms because of the cumbersome terms of Taylor series expansions. In 
comparison to other approaches, our formulas are rather transparent and forward, 
without losing the precision apparently. The mathematical tools are explained briefly 
in section 2, while the formulas are given afterwards. In section 3 and 4 we intend 
to inspire undergraduate students by some typical examples mentioned above after the 
derivation of the approximate formulas. Section 5 is devoted to our conclusions. 

2. Mathematical preliminary: methods and formulas 

The Lagrange inversion theorem [4] is a remarkable tool famous for its ability to give 
explicit formulas where other approaches run into stone walls. Lagrange-Biirmann 
formula [5], a special case of Lagrange inversion theorem, has found many applications, 
such as evaluating roots of certain transcendental equations and obtaining expansions 
of a function in powers of a related but different function. 

Let the function f{z) be analytic in some neighborhood of the original point z = 0 
of the complex plane with /(O) ^ 0 and satisfy the equation [4] 
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There exists two positive numbers a and b such that for |t(;| < a the equation(3) has 
just one solution in the domain \z\ < b. Lagrange-Burmann formula tells us that the 
unique solution is an analytical function of w with [5] 


z = 


E 

n=l 


r 

n! \-dz'^~^ 



(4) 


In contract to Taylor expansion, Bade approximat [6] has abundant applications 
in physics because of its fast convergence speed and elegant form. Bade approximant 
is a type of rational fraction approximation to the value of a function. This structure 
of approximant enables an effective reconstruction of function’s singularities over the 
whole range using its series expansion obtained for small values of its variable. The 
\p, q] Bade approximant denotes a fraction of polynomial Pp{x) = of degree at 

most p, and polynomial Qq{x) = ^nx"' of degree at most q 

The fraction consisting of polynomials Pp{x) and Qq{x) which approximates function 
f{x) = Yhn fnX"' is determined by the equation [6] 

m - ^ ( 6 ) 

where the symbol 0{x^) stands for the value of the order x”. go = 1 is always assumed 
for convenience. It is known that in many cases a higher accuracy of approximation 
is achieved for fraction of polynomials of identical degree, thus the degrees of both 
numerator and denominator are set to be 2 hereafter. 


3. Selected applications of transcendental equation (1) 


3.1. Fade approximant to (1) 


We plan to hunt for approximate formulas relating to (1) by the Lagrange inversion 
theorem and Bade approximant in the current section, while the approximate formula 
of (2) will be provided in the next section. 

The analytical solutions to (1) are absent until now. From the perspective of 
graphical method, the solution to (la) (or (16)) is equivalent to the solution of the pair 
of equations y = kx and tangent function y = tanx(or cotangent function y = cotx). 
The discussion will be limited to non-negative solutions because physical quantities 
involved in current paper only relate to positive roots. Suppose that the zero in the 
domain [0, ^tt) is denoted by Xq while the zero in the domain [{n — 4)7r, (n -|- \)ti) is 
denoted by x^ with integer n > 1 for (la) and (16) respectively. 

Let us concentrate on the hrst zero of (1) hrstly. The hrst zero of (la) is trivial(i.e. 
x^ = 0) for any k, while situation is far more complicated for (16). The [2,2] Bade 
approximant yields [15] 


X 


0 


1 1 -L I 103 ,^-2 

1 ^ 4044 '^ ^ 5593 


K 1 -I- 655 -1 I 255 -2 

^ 1348 ^ 3704 


1 ^2 +0.3192k+ 0.0184 
k2 + 0.4859k + 0.0688 






(7) 
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for big K, while for small k it is 


TT 1 + 2ft; ~\~ 

° ^ 2l + 3K+(2 + fJ)K2' 

(8) 

We pass now to consider the other roots. By means 

of Lagrange-Biirmann formula 

(4), we arrive at (see appendix A for detail) 



(9a) 

~ /3n0+(l//5n)) 

(96) 

where «„ = (u + l/2)7r, fdn = nn, and 0±(x) are defined 

as 

4>±{x) = 1 ± -x^ - ^x'^ + 0{x^). 

(10) 

Since (j)±{x) are even functions, the [2,2] Fade approximant to (10) turns out to be 

, , , (6ft; ± l)x^ ± 3ft;^ 

(11) 

Specially, if ft; = 1, the solutions to (1) are 


, 3q!^ — 5 

T' rv ''' 

^ ”3a2_2’ 

(12a) 

- n A ^ 

~^"3/32+4‘ 

(126) 


We mention here that Frankel [16] once obtained a fairly accurate approximation 
solution to (la) in terms of the Taylor series expansion of arctanx. He reached at [16] 


xj) ~ Q!„ — (1 H—-)arccotaji. (13) 

We note that (13) is sightly different from the original version derived by Frankel himself 
because of the fact that arctanx + arccotx = 7r/2 if x > 0. Motivated by his inspiring 
thoughts, we manage to hnd the counterpart formula for {lb), namely 

1 + 

Xn - f3n+ ^ arccot/ln. (14) 

' Pn 

Detailed comparisons among our formulas and Frankel’s formula will be made to 
end this subsection. Table 1 and table 2 show the accurate values x^ and the absolute 
errors calculated by (9), (12), (13) and (14) with n = 1. It can be distinguished that the 
formulas based on Fade approximant have an advantage over Frankel’s formula sightly, 
but the advantage will be extended with the increase of the order of Fade approximant. 
Besides, both tables tell us that Fade approximant is around twice as precise as Taylor 
series expansion, which, of course, suggests the superiority of Fade approxiant. 
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Table 1. We present the first 10 non-trivial roots of (la) with k = 1. It can be 
distinguished from the column of x^/tt that x+ approaches to (n + l/2)7r as n —>■ oo. 
The errors of three different formulas, which are in the unit of 10“^, are shown in the 
last three columns in order. _ 


n 

Exact values 

X + 

Errors(xl0 ^) 

TT 

Fade 

Frankel 

Taylor 

1 

4.49340946 

1.43029665 

0.20508427 

0.45855420 

0.40225822 

2 

7.72525184 

2.45902403 

0.01474265 

0.03420796 

0.02977709 

3 

10.90412166 

3.47088972 

0.00268424 

0.00629142 

0.00546478 

4 

14.06619391 

4.47740858 

0.00075749 

0.00178279 

0.00154718 

5 

17.22075527 

5.48153665 

0.00027654 

0.00065221 

0.00056576 

6 

20.37130296 

6.48438713 

0.00011965 

0.00028254 

0.00024503 

7 

23.51945250 

7.48647425 

0.00005841 

0.00013804 

0.00011969 

8 

26.66605426 

8.48806870 

0.00003121 

0.00007378 

0.00006397 

9 

29.81159879 

9.48932662 

0.00001788 

0.00004229 

0.00003667 

10 

32.95638904 

10.49034444 

0.00001084 

0.00002564 

0.00002222 


Table 2. We present the first 10 non-trivial roots (neglecting Xq) of (16) with k = 1. It 
can be distinguished from the column of x~/nir that x~ approaches to nir as n —>■ oo. 
The errors of three different formulas, which are in the unit of 10“^, are shown in the 
last three columns in order. _ 


n 

Exact values 


Errors(xl0 ^) 

riTv 

Fade 

Frankel 

Taylor 

1 

3.42561846 

1.09040822 

-0.36000169 

-0.18196111 

-0.87179656 

2 

6.43729818 

1.02452783 

-0.01575732 

-0.00868217 

-0.03331845 

3 

9.52933441 

1.01109378 

-0.00222643 

-0.00124921 

-0.00458176 

4 

12.64528722 

1.00627998 

-0.00054192 

-0.00030606 

-0.00110449 

5 

15.77128487 

1.00403118 

-0.00017970 

-0.00010180 

-0.00036460 

6 

18.90240996 

1.00280399 

-0.00007269 

-0.00004125 

-0.00014712 

7 

22.03649673 

1.00206211 

-0.00003376 

-0.00001918 

-0.00006823 

8 

25.17244633 

1.00157982 

-0.00001736 

-0.00000987 

-0.00003505 

9 

28.30964285 

1.00124880 

-0.00000965 

-0.00000549 

-0.00001947 

10 

31.44771464 

1.00101185 

-0.00000571 

-0.00000325 

-0.00001151 


3.2. Selected applications of transcendental equation (1) 

Effect of spring mass on the frequency of oscillator The oscillation of a spring-mass 
system, where a spring is suspended vertically and a mass m is hung from the bottom 
end of the spring, is a rather classical topic commonly studied theoretically and 
experimentally in introductory physics courses. The consideration of the correlation to 
the spring oscillation due to the non-ignorable mass of spring has led to many papers, 
and the empirical law that 1/3 the mass of the spring should be added to the mass of 
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the hanging object is frequently-quoted(see reference [7] and references therein). The 
motion of this celebrated system is governed by [17] 


{ ^tt ® ^XX 0 ) 

u\x=o = + Y^Utt)\x=i = 0, (15) 

u\t=o = ^X,Ut\t=0 = 0 , 


where Y, s and p are the Young’s modulus, cross-sectional area and mass density of the 
spring, and a = \/Y/p is the wave speed in the spring. The solution to (15) is 


CX) 

u{x,t) = 

n=l 


2ps 


]\ —^^0 cos(y^at) sin(yA^: 


mps 


(16) 


where satishes the transcendental equation 

cot(\/A/) = — y/Xl. (17) 

mo 

Taking into account the Hook’s law in the language of Young’s modulus, which is 
dehned as the ratio of the stress(force per unit area) along an axis to the strain(ratio 
of deformation over initial length) along that axis in the range of stress, we obtain the 
relation Y = kl/s where k is the stiffness of ideal spring, and therewith (17) becomes 


cot 




(18) 


where the angular frequency uj = a\/A. After a trivial trick, the frequency of the spring 
turns out to be 


where 


u = 


/ ^ 

\l m + ^mo 

(19) 

1 


~ r?2 

it(r^) = rp 

(20) 


with r = m/mQ. It can be concluded from (16) that the motion of the spring is 
the superposition of inhnite numbers of simple harmonic vibrations. However, the 
amplitudes of all the other vibrations are much smaller than the principal vibration’s. 
Many authors, such as Christensen [7], showed explicitly that the small r case necessarily 
transitions to the large r case so that the lowest normal mode is always the dominant 
one. In this occasion, the fundamental frequency of the spring can be calculated by (7) 
or (8) subsequently. 

Now let us consider two extreme cases. If r is large enough, then rj :^(1 ~ 
therefore ^min = lim r[(l-A)"2_l] = 1 On the contrary, = hm[4(l-r)"^ _r] = 

1 -l-CXD 1 -l-O 

since ?7 ~ |(1 — r) when r is a small number. Specially, if the mass of the spring 
and hung object are equal, namely r = 1, one obtain p 0.86 from (7) and therefore 
hnd the effective mass coefficient is ^ ~ 0.35, which is in accordance with the fact that 
I Y ^ Y ^ and satishes the empirical law mentioned above. 
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Fraunhofer single slit diffraction One of the most long-standing problems of classical 
optics is the various types of diffraction. The theory of Fraunhofer single slit diffraction 
predicts that the spatial pattern of light intensity on the viewing screen by a light wave 
passing through a single rectangular-shaped slit is given by [8] 


I = Ii 


sin^ u 




( 21 ) 


where u = 7r6sin6*/A and Jq is the light intensity at 6* = 0°. The hrst and second 
derivative of I with respect to M(or 6 more precisely) are badly needed in order to hnd 
the maxima and minima of the diffraction pattern. The minima occur when u = nvr, 
n = ±1, ±2, • • •, while the maxima exist on condition that 


tanu = u. (22) 

The trivial zero of (22), i.e. u = 0 corresponds to the primary maximum of the diffraction 
pattern, while the non-trivial zeros indicate the secondary maxima. We see that the 
secondary maxima are not exactly half way between any two adjacent minima, they 
occur slightly earlier and move toward the center with increasing shift u. Figure below 
shows the variation of the intensity distribution with the distance. 



Figure 1. The diffraction pattern of Fraunhofer single slit diffraction. It can be 
distinguished that intensities of secondary maxima are dramatically less than the 
principal maximum, and the positions of secondary maxima are not exactly half way 
between the adjacent minima. 

Furthermore, intensities of these secondary maxima are much less than primary 
maximum and fall off rapidly as moving outwards. The relative intensity of the n-th 
secondary maximum to the primary maximum can be determined from (12 a) as 

f 1 q 1 

f = o(^rT^--) (23) 

/o 2 3 q!^ — 2 

where «„ = (n -|- l/2)7r. For example, the relative intensity of the hrst three secondary 
maxima are only 4.7%, 1.7% and 0.8% of the principal maximum, respectively. 
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Schrodinger equation with single 6 potential We shall discuss a problem encountered 


in many introduction course of quantum mechanics in detail. Now let’s consider a 
microscopic particle confined in an 1 -dimensional infinite square well in 0 < x < a, and 
assume that a 5-function interaction is added in the middle of the interval with constant 
strength 7 , whose sign indicates the interaction is attractive or repulsive. Therefore, 
the potential of the particle is [10,15] 



0 < x < a, 
otherwise 


(24) 


The wavefunctions can be classified according to their symmetries under the 
transformation 'ijjn{x) = ±'0n(fl — x) since the Hamiltonian is symmetric with respect to 
the transformation x ^ a — x. Thereby only the even parity eigen-states ipnix), which 
satisfy the Schrodinger equation 



(25) 


are influenced by the additional interaction. At point x = a/2 the wavefunction should 
be continuous but its derivative makes a jump proportional to the strength of the 6- 
function interaction, i.e. 



(26) 


with e —)■ 0. In general, the energy En relates to the quasi wave-vector /c„, that satishes 
the energy eigenvalues condition: 



tau (-) = --(-) (27) 


where b = 2fF/m'y is the characteristic length of 5-potential. The energy are positive 


if the interaction is repulsive, situation will be more complicated if the interaction is 
attractive. In the case of 7 < 0, ground state energy for instance, will go through a 
process from positive to negative with increasing intensity of the interaction. However, 
unlike the other eigenstates, the E = 0 energy eigenstate of the Schrodinger equation is 
not a sinusoidal function. Instead, the time-independent Schrodinger equation simplihes 
to d’^'ijj/dx'^ = 0 inside the well, and yields ^|J = Ax+B, where A and B are constants [15]. 
The energy eigenstate has the proper discontinuity in its slope at the middle point of 
the well such that [ 10 ] 




0 < x < a/ 2 , 


(28) 


The combination of (28) and (26) yields the critical intensity 70 
(27) is reduced to 



and therefore 



( 29 ) 
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Equation (12a), readily, tells us that the solution of (29) satisies the identity {kna/2) — 
(u7r/2)^ “2(1 + 2/3(n7r)^) with n = 3,5, 7, ••• approximately, then the even-parity 

energy levels of the system considered are 


El = 


HVn 

2m 


= E® + 

^ 4a' 


n = 1, 


Ef+ 2(1 + —n = 3,5,7, 
V 3(n7r)"'/ a 


(30) 


2 2 fe 2 

where En energy of one-dimensional inhnity square well without <5- 

potential. The errors of the energy obtained from (30), which can be inferred from 
table 1, are rather small. Furthermore, (30) suggests us that the 5-potential has a 
more signihcant influence on ground state energy than the excited state energy since the 
prefactor 7r^/4(!^ 2.5) for n = 1 is larger than 2(l-|-2/3(n7r)^) ~ 2 that for n = 3, 5, 7, • • •. 


4. Selected applications of transcendental equation (2) 

4 . 1 . Fade approximant to (2) 

The Lambert W function is a multivalued function with complex variable though, of 
special relevance to scientihc applications are the solutions when the argument is purely 
real number. The two real solutions of (2) are the branches ITo and lT_i where Wq is 
the principal branch of the W function. 



X 


Figure 2. Two real branches of the Lambert W function. Red line: Wo{x) called 
the principal branch is defined for — 1/e < x < -l-oo. Blue line: W-i(x) is defined for 
— 1/e < a; < 0. The two branches meet at the green point (—1/e, —1) [18]. 


The principal branch Wq is analytical at a: = 0 and its Taylor series expansion [13,14] 
^ - x"' = x(^l - X + -b + 0(x^)^ (31) 

n=l 


can be obtained in light of Lagrange inversion theorem. The [2, 2] Fade approximant 
reads 



= X 


1 I 19™ I 17 ™2 

W -r ggX 

1 _L I iOl.r.2 

1 -t- qq X 


+ C>(x=). 


(32) 
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After rounding the fractional coefficients by the special method introduced in [6], (32) 
is reduced to 


-nrO)/ N 3 + 6x + x‘^ 

Wq (x) = ^ ^ + 0 {x^). 


(33) 


3 + Ox + 5a;2 

which, of course, is not only more elegant but isn’t accomplished at the sacrihce of 
accuracy. In fact, the substitution x —)■ ln(l + a;) of the pre-factor is a good succedaneum 
since the slope of W{x) around zero is no more than 1. We therefore arrive at 


Wq^^\x) = ln(l -I- x) 


1 + + 21^2 


10 " 


1 _|_ 143 rp I 713 rp2 

-L ~ Af\ vt/ \ 


+ 0(1= 


(34) 


40 ' 240" 

if we utilize the fact that ln(l-|-a;) = Readers can refer to the appendix 

below for more details. Similar method impels result such that 


T77(II)/ N , .2 + Qx + 4:X^ ^ 

Wq (x) = ln(l a^)7yTWrTW3 + 


2 -h 7x + 6x'^ ' \ / ( ) 

Figure 3 shows the precision of the two formulas, i.e. (32) and (34) introduced 
above. It can be concluded that Fade formula of type II is superior to type I 
overwhelmingly, especially for positive argument. Therefore, we recommend that (34) 
is the perfect approximate formula. 



Figure 3. The curves of the logarithm of the absolute errors(A = Ig (| I)) 

regarding to the two different formulas in the interval —l/e<x<l are plotted where 
W^(x) and Wq{x) are the approximate and exact values of the W function. The 
dashed line and the solid line, from the top down, represent Fade formula of type I 
and II respectively. 

Last but not the least, let us focus on an equation relating to Lambert W function. 
The equation is 

= a(x — b) (36) 

where a, b, c are constants. A great deal of quantities, such as the decay constant of 
an exponentially decaying process [19] and the time constant of a process subject to a 
linear resistive force [11], satisfy the same equation as shown above. The solution to 
(36) is 

x = b + -W(-e-‘^’^). 
c a 


(37) 
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Schrodinger equation with double 6 potential Let’s extend to the case of a particle in 
a double delta function well which can be used to describe the behaviour of electronic 
terms of Mathematically, the potential is described by [9] 


V{x) = + a) + S{x — a)j 


(38) 


where 7 , a > 0. We wouldn’t repeat the tedious derivation again, instead, we recommend 
the readers to some professional literatures, such as [9] for details. Here we just sketch 
the general conclusions. The potential is an even function, so all the solutions can be 
expressed as a linear combination of even and odd solutions, namely the so-called even 
parity and odd parity. The bound energy of the particle, in either situation, can be 
expressed as 

2 m 

where the wavevector is the solution to equation 

= (40) 



with b = 2 ^, and the plus and minus sign correspond to the even parity and odd parity 
solution, respectively. Thanks to (37) we manage to obtain the bound state energy 


^ ~8ma^ ib^ ^^b^ 


(41) 


Therefore, it can be inferred that the bound state energy always exists for even parity, 
while for odd parity it only exists on condition that a > b. In other words, there’re 


two bound states if the intensity of the interaction 7 > 2 ^, otherwise there exists only 
one bound state. Our approximation formula (34) is sufficient to determine the energy 
since the argument of the Lambert W function in (41) is bounded from —1/e to 1/e 
for both cases. For example, if a = 2b, the bound state energy are —0.6148^^ and 
-0.3176-^ respectively, while the accuracy values presented in [9] are -0.614782-^ 
and —0.317454^^ respectively. The fact that the relative errors are 0.003% and 0.05% 
or so indicates that our approximate results are in glorious agreement with the exact 
ones. The exchange energy IS.E dehned as the difference between and E~ can be 
determined subsequently. 


Wien’s displacement law Planck’s seminal work for blackbody radiation, whose result 
conflicted dramatically with classical mechanics, opened the era of quantum theory. The 
concept of blackbody radiation, along with the associated Stefan-Boltzmann law and 
Wien’s displacement law, is a crucial pillar of modern physics. The Planck spectral 
distribution is given by [19] 


P(A,T) 


87rhc 1 

^5 g/ic/A/cT _ 


(42) 
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where A is the wavelength, T is the temperature, c is the speed of light, and h and k 
are the Planck and Boltzmann constants respectively. Wien’s displacement law gives 
the wavelength at which (42) has the maximal intensity. To hnd the extremum of the 
specihc intensity, it’s necessary for us to take the derivative of (42) with respect to 
wavelength A and utilize the often-used unitless variable x = hc/\kT. Apparently, we 
have the fancy equation 


(5 - x)e^ = 5 


(43) 


whose non-trivial zero is of great interest. In order to solve (43) analytically, several 
attempts have been made to express the root in integral representation and series 
representation over the past one century. Siewert and Burniston found that the 
analytical solution of (2) relates to the canonical solution of the Riemann problem [20], 
and Siewert applied it to Wien’s displacement law to observe that [21] 


xo = 4 exp 


f 1 

POO 

V 71 j 



arctan 


TT 


TT 


dt 

t + 5 


(44) 


dn 5 — 5 — t — Int/ 

Luck and Stevens, at the meantime, presented another integral representation by virtual 
of Cauchy’s integral theorem and some basic concepts of complex integration [22], i.e. 


w{6)e^^^d6 

Xn = 5-^4 - 


(45) 


where w{9) = i other hand, Andersen [23] and Vial [19] suggested 

that only the hrst three terms of the Taylor series expansion in (31) can reach a 
satisfactory precision. 

As a matter of fact, (43) is nothing but (36) with a = 1/5, 5 = 5 and c = 1. 
Therefore Wien’s displacement law can be obtained with an elegant expression for 
displacement constant b = hc/kxQ where 


Xo — 5 + Wn(-de ^). 


(46) 


It can be distinguished from (44), (45) and (46) that the nontrivial zero expressed 
by Lambert W function pushes the mathematical structure of the law into the 
most comfortable territory, and the calculation based on it is especially accuracy. 
Furthermore, our approximation formula (34) gives Xo = 4.965114231797, which owns 
ten decimal places precision compared to the accuracy value Xq = 4.965114231744. The 
relative error, without doubt, is rather small, and it is of the order 10“^^. What’s more, 
even the simplihed version (35) can uncover four decimal places precision. 


5. Conclusion 

In spite of the fact that the solutions to many transcendental equations, such as (1) and 
(2), can be expressed analytically in power series or named special functions, but it’s 
not really convenience for instructional purpose in classroom. This paper has provided 
a detailed illustration of how Lagrange inversion theorem and Bade approximant can be 
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applied to solve transcendental eqnations analytically, and express the root(s) of those 
eqnations in the form of ratio of rational polynomials. Furthermore, our approximation 
method should be quite powerful in handling with the non-transparent aspect of some 
rigourous results, which could eventually be expressed in more brief forms, without 
losing the physical interpretations. Traditional graphical method, however, is not only 
troublesome but too oversimplihed to reach desired accuracy. While all of the drawbacks 
can be overcome by the built-in algorithms in the softwares such as Matlab and 
Mathematica, the mathematical softwares tend to give numerical results, or just make 
the so-called closed-form solutions be lengthy, which are useless in practice. Our method, 
therefore, is superior in teaching because it is easy to handle and owns satisfactory 
accuracy. An excellent agreements were achieved between our approximate results and 
analytical formulas, which prove that those formulas are rather precise in practice, and 
therefore can act as valuable ingredients in the standard teaching curriculum. 
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Appendix A 


we begin the process of putting the equation(la) into the form required by the Lagrange 
inversion formula hrstly. The asymptotic expression ~ {n + ^)'k holds as n ^ oo. 
Let us set = w~^ — x where w = then (la) is reduced to 

tan(x) = tan -|- -)7r — = cot( 2 ;) = — z). (A.l) 

Therefore, (A.l) can be arranged to (3) with 


/+(^) = 


2 ;(cos z + nz sin z) 


(A.2) 


Ksmz 

Similarly, we can get the expression for /_. We thus can obtain the approximate formulas 
to (1) according to (4) since f±{z) is analytical around z = 0 and /(O) = 1/k. 


Appendix B 


The calculation of W^^\x)^ where the pre-factor is ln(l Fx) rather than x, is not really 
straightforward. To begin with, let us dehne the auxiliary function 


M{x) 


W{x) 
ln(l -I- x) 


n=0 


(B.l) 
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where an are the Taylor expansion coefficients, if we ntilize the fact that ln(l + x) = 
and the eqnation (31), we hnd that 


1 11 2 43 


2651 


M{x) = 1 — -X + —X — —x'^ + X + 0{x^) 


(B,2) 


2 " ■ 12 " 24 720 

after rearranging the coefficients of the same power. We therefore can get the second 
order Fade approximant of the Lambert W fnnction aronnd zero as is shown in (34) by 
virtne of the Fade approximant of (B.2) 


M{x) = 


1 _L 1^^ I ^.^2 

+ +0{x^). 


1 _1_ IMt- 4- 

W 40 W 240-^ 


(B.3) 
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